g-Breathers in Discrete Nonlinear Schrodinger arrays with weak disorder 



OO 

o 
o 

(N 
o 

Q 



> 
in 

(N 
OO 

o 



M. V. Ivanchenko 

Department of Applied Mathematics, University of Leeds, LS2 9JT, Leeds, United Kingdom 

Nonlinearity and disorder are key players in vibrational lattice dynamics, responsible for localiza- 
tion and delocalization phenomena. g-Breathers - periodic orbits in nonlinear lattices, exponentially 
localized in the reciprocal linear mode space - is a fundamental class of nonlinear oscillatory modes, 
currently found in disorder-free systems. In this paper we generalize the concept of g-breathers 
to the case of weak disorder, taking the Discrete Nonlinear Schrodinger chain as an example. We 
show that g-breathers retain exponential localization near the central mode, provided that disorder 
is sufficiently small. We analyze statistical properties of the instability threshold and uncover its 
sensitive dependence on a particular realization. Remarkably, the threshold can be intentionally 
increased or decreased by specifically arranged inhomogeneities. This effect allows us to formulate 
an approach to controlling the energy flow between the modes. The relevance to other model arrays 
and experiments with miniature mechanical lattices, light and matter waves propagation in optical 
potentials is discussed. 

PACS numbers: 6.3.20.Pw, 63.20.Ry, 05.45.-a 



A wealth of physical systems from natural crystals 
to the cutting-edge technology products like micro and 
nanomechanical system arrays are spatially extended and 
discrete. Interaction between their elements is a key 
source for a number of fundamental dynamical and statis- 
tical physical phenomena including thermal conductivity, 
wave excitation and propagation, electron and phonon 
scattering. To provide with a full understanding of these 
processes the theory of collective vibrational modes is 
in demand. The principal question to be answered is 
the effect of the two fundamental features of any lattice: 
nonlinearity and disorder. 

In recent decades we have witnessed a remarkable 
progress in studying their impacts separately. Nonlinear- 
ity induces interaction between linear normal modes and 
energy sharing if it is strong enough (the Fermi-Pasta- 
Ulam (FPU) problem) or exponential localization 

of exact periodic solutions (discrete breathers) in space 
01 • Independently, disorder leads to exponentiallylocal- 
ized linear vibrational modes (Anderson modes) 

However, the concurrent effect of nonlinearity and 
disorder has not received a satisfactory full description 
yet. Strongly disordered and weakly nonlinear lattices 
enjoy intensive research, in particular, on continuation 
of Anderson modes into nonlinear regime wavepacket 
spreading Ifl], light propagation in photonic lattices [?], 
and Bose-Einstein condensate (BEC) localization in ran- 
dom optical potentials [8|. In contrast, little is known in 
case of pronounced nonlinearity and weak disorder. Im- 
portantly, this situation is reaUzed in micro and nano- 
mechanical oscillatory arrays that are often driven into 
nonlinear regime, while the spatial disorder is constantly 
reduced by improving fabrication techniques [^, [lo| . 
On the atomic scale, the surface vibrational modes are 
thought to be a source of selective catalytic properties 
of three-dimensional gold nano-clusters for a variety of 



namics in random optical media are equally strong mo- 
tivating problems. 

g-Breathers (QBs) present a recently discovered funda- 
mental class of nonhnear oscillatory modes. They are ex- 
act time-periodic solutions to nonhnear lattice equations, 
continued from linear normal modes and exponentially 
localized in the linear mode space. Introduced to explain 
the FPU paradox (energy locking in the low-frequency 
part of the spectrum, recurrencies, and size-dependent 
stochasticity thresholds) they have been found in 

two and three dimensional FPU arrays and discrete non- 
linear Schrodinger (DNLS) lattices last year quan- 
tum QBs were observed in the Bose-Hubbard chain [l^ . 
QBs have been suggested as key actors in a BEC pulsat- 
ing instability jl5| a nd a four- wave mixing process in a 
nonlinear crystal 16| . 

In this paper we extend the concept of g-breathers to 
random arrays, implementing the paradigmatic DNLS 
model as an example. The cornerstones of our approach 
are continuation of QBs into non-zero 'frozen' disor- 
der, taking a nonlinear locahzed solution as a seed, and 
analysing statistics then. We show that QBs display the 
crossover from the exponential localization near the cen- 
tral mode to the power-law decay at a distance. Their av- 
erage linear stability threshold in nonhnearity keeps the 
same value in the first order approximation. The vari- 
ance increases linearly with disorder, manifesting high 
sensitivity on particular realizations. Finally, we demon- 
strate, that the superimposed periodic modulation of the 
linear coupling strength can be a means of the energy 
flow control. 

The DNLS lattice is represented by the Hamiltonian 

H = + DKn)^Pn-ltpn + (1 + D Kn+l)tpn+l'tpn + 
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and the equations of motion are i^pn = dH/diJj^: 

iljjn = {l+DKn)ljjn-l + ii + DKn+l)lpn+l+fi\lpn\'^ Ipn (2) 

Here V is a complex scalar that may describe the slow 
small-amplitude dynamics of a classical nonlinear oscil- 
lators array 171, probabilit y am plitude of an atomic 
cloud on an optical lattice site [19||, or the amplitudes of 
a propagating electromagnetic wave in an optical waveg- 
uide |20|. Zero boundary conditions apply: ipo = i1^n+i = 
0. /i and D are the nonlinearity and disorder parameters, 
random k„ G [—1/2,1/2] are uniformly distributed and 
uncorrelated: (K„Km) — (r^Sn^m, = 1/12. Beside the 
total energy, the norm B — Y^ IV'nP is conserved. Chang- 



ing /i is strictly equivalent to changing the norm B, thus 
we fix i? = 1 further on. The canonical transformation to 
the reciprocal space of normal modes with new variables 
Qq(t) is given by 
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FIG. 1: The average mode energy distribution in QBs with 
increase of disorder, where = 0.1, N = 100 : (a) the low 
frequency mode qo = 11 and (b) the middle frequency mode 
go = 53. Filled circles are analytical estimates l|lip 



and the dynamics in this space reads: 

p,r,s 



where p 
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(4) 
are 
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the normal mode frequencies for the linear disorder 
free system with fj, — D — 0. The nonlin- 
ear intermode coupling coefficients are Gq^p^r.s - 



^(_l)(±p)(±.)(±.) {Sq±p±r±s,0 
± 

the disorder induced ones read Vr 



q±p±r±s,±2{N+l) 
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N-l 



Tvq{n+1) ^. TTpn \ 
^^"^ N+1 ^^"^N+l' 



E «:„(sin^sin2£^ 

n—l 

In the disorder-free case QBs are time-periodic sta- 
tionary solutions 'ipn{t) = 4>n exp(irit) with the frequency 
Q. and time-independent amplitudes (/)„ localized in nor- 
mal mode space. In the g-space they have the form 
Qq{t) = Aqex\){int), the amplitudes of the modes Aq 
being time-independent and related to the real-space am- 
plitudes by the canonical transformation ([3]); the mode 
energies are defined as Bq — \Aq'^ . Here we focus on 
time-reversible periodic orbits and, thus, consider Aq to 
be real numbers. The amplitudes satisfy a closed system 
of algebraic equations: 



(f2g — n)Aq — — ^ ^ Gq^p^r,sApArA^ + d ^ ^ Vq^pAp, 



p,r,s 



J2\Aq\'~B^0 



(5) 

Our methodology consists of two steps. Firstly, we 
take a known QB solution for non-zero nonlinearity [l3 |. 
A particular realization of {k„} is chosen and d regarded 



as the disorder parameter. Together with the nonlin- 
earity parameter p, it is assumed to be small p,d 1. 
Then, an asymptotic expansion in powers of {p, d} is de- 
veloped. Subsequent Hnear stability analysis employs the 
constructed solution. Secondly, statistical properties of 
the QB solution and the instability threshold are ana- 
lyzed. 

Continuation of QBs from ^ 7^ 0, = to ^, Z? 7^ 
exploits the same ideas as from p — D — Oto^^ 
Q,D — Q For small amplitude excitations the non- 
linear and disorder terms in (U can be neglected and 
the g-oscillators get decoupled, their harmonic energy 
Bq — \Qqf' being conserved in time. Single g-oscillator 
excitations {Bq ^ for q ^ qo only) are trivial stationary 
and (^-localized solutions for (3 = D = 0. 

In the disorder-free case such periodic orbits can be 
continued into the nonlinear case at fixed total energy 
[l3 | by solving the system of algebraic equations |(5]), 
granted by the implicit function theorem [2l|, as the non- 



resonance condition flq^ ^ ^q=iqo holds. This is vaHd for 
d <C 1 as well, for the spectrum remains non-resonant 
with the probability 1 Numerically, we were able to 
continue QBs into the (3,D ^ domain for all parameters 
taken. 

Typical results for the low-frequency and middle- 
frequency QBs are shown in Fig[TJ They demonstrate 
the crossover between the exponential localization and 
the disorder induced background. The disorder-free ex- 
ponential localization persists in some neighborhood of 
the central mode for sufficiently small disorder, but is 
range shrinks as disorder grows. High-frequency QBs be- 
have analogously. 

Let us construct an asymptotic expansion for the 
QB solution. We assume p, d <C 1 and start from 
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the disorder-free QB profile A^^ for the modes q 
3qQ,. . . ,(2n + 1)90)- . ■'ti N in the leading order of p [li 

/i(7V + 1) 
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We seek an asymptotic expansion in powers of c? ^ 1: 
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0{(P,pd), where A^"' = A^ 



tion into |[5|) gives 
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The ensemble average of the "disorder contribution" to 
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the energy 
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that approximates well the disorder-dominated part of 
the numerically obtained QB profiles in different parts 
of the linear spectrum (Figd]). The power-law decay 
{B^'^) oc {q ~ qo)~^ fits in the large part of the q- 
space. One can estimate the crossover location between 
the exponential decay and the power-law, in particular, 
when the modes next the to the central one become ex- 
cited almost equally well. Letting q = qo + 1 we obtain 
the "small" a^D < n/{2{N + if/^) and the "large" 
a^D ^ tt/ \/ {N + 1) disorder criteria for delocalization 
of the least robust modes qo = 1,N and the most ro- 
bust one qo = N/2. Thus, single-cite centered modes do 
not exist above the size-dependent threshold in disorder 
magnitude. (Note, that delocalization in the mode space 
approximately corresponds the onset of the Anderson lo- 
calization in the direct space.) 

The linear stability of QBs is determined by con- 
sidering the evolution of small complex-valued pertur- 
bations (q{t) to the stationary solution Qq{t) = 
{Aq + Cq{t)) exp (iflt). In linearized equations the stabil- 
ity requires all the eigenvalues be negative. Numerically 
we solve the correspondent problem in the direct space 
In the following we restrict our attention to the 
low and middle-frequency QBs, leaving the more com- 
plex case of qo > N/2 (when for D = the instability 
threshold behaves erratically vs. 90 for the future 
study. 

We find, that the instability develops similarly for zero 
and non-zero disorder, the increase or decrease of the 
threshold p* sensitively depending on a particular real- 
ization (Figl2]). The average (/i*) remains very close to 
the zero-disorder value Po- In contrast, the variance a^- 
is significantly growing, depending on D almost linearly 
(FiglSl deviations being observed when the probability of 
/i* being next to zero becomes substantial). 
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FIG. 2: (a) The maximal eigenvalues 9 of QBs with qo — 
15, A'^ — 128 and two different sets of {«;„} vs. the nonlinearity 
coefficient for several values of disorder strength D. For one 
realization of disorder the instability threshold is increasing 
with D, for another - decreasing, (b) The variance of the 
QB instability threshold a^* vs. D. Solid lines are analytical 
estimates l|lip 



The analytic study of the QB stability has not been 
done before (even for D = 0) and we present it here 
for the first time (restricting to go < N/2 as above). 
Linearized equations for small perturbations read: 



iQ — (f2 — flq)Cq + ^Bq„ ^ Gq^qg^qg^p{(* + 2(p) 



(9) 



In analogy to the FPU chain we suggest (and verify 
that by comparison with the numerical results) that the 
eigenvectors for the main instability be almost parallel to 
the subspace {Q = : q ^ qo±l}. Thus we arrive at a 
simpler task of finding eigenvalues of the system of two 
complex-valued linear equations (retaining 0{p, d) terms 
only): 



iCqo-l = (f^ - ^qo-l + dVqo-l,qa-l)Cqo-l + 

+ ^PBqaiCqo + l + 2Cgo + l) + ^Vg,, + i Cqo +1 , 

*C<jo + l = (f^ - ^qo + 1 + c'K?o + l,go + l)Cgo + l + 
1 

+ 2 



(10) 



pBqoiCqo-l + 2C90-1) + c'^^Q-l.go + lCgo-l 



After an extensive algebra one finally gets the bifurcation 
point: 
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where AK,, 



90-1,90- 



2Kjo,9o + ^0+1,90+1; and 
It shows a good coin- 



the disorder-free /Zq = 2_b (jv+i) • 
cidence with the numerical results (Figl2]). Note, that in- 
creasing the chain length decreases the instability thresh- 
old and increases its variation. Thus, in sufficiently large 
arrays the solution will loose stability at very small non- 
linearities with the probability, almost equal to that of 
AVqg^qg being negative, which is 0.5 in our case. 

If the instability depends that sensitively on the 
disorder realization, there must be certain classes of 
inhomogeneities that augment it or suppress. Iden- 
tifying them offers the possibility of controlling the 
energy flow in the mode space by designing specific 
impurities and, further, changing them in time. The 
disorder induced correction in (fTTj) reduces to AVq^^qg — 



7rgo(2w+l) • 2 7r(2Tt+l) 
COS sm 2(N+1) 



0{N- 



N+ 

Note, that it is Hnear in k„, and, therefore, one can 
represent k„ as a sum of spatial Fourier components, 
their contributions being additive to A V^^ . 

Let us consider a harmonic inhomogeneity k„ = 
^ cos ( '^^^j^^y +^), where ip is the phase shift. It is 
natural to expect the absolute extrema of A.Vqg^qg{p,ip) 
(and the maximal gain or loss in stability) to be reached 
for p = 2qo. This case yields AVqg^qg « VN + Icost^. 
Thus, the bifurcation point reaches its maximum and 
minimum for ip = ir and (p — respectively, giving 

D{N+1) 

zero shift for ip = zL^. Analogously, for p = qo zL I one 

8-D(jV+l) 
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At the same time one gets a 



gets ^ 



1 ± 



■ sm(p 



which is — f shifted 



in p and has a bit smaller amplitude. For p ~ qo zt2 it 
reads ^* ~ (Iq ( 1 + cos p ) , which is tt shifted in 

\ 27r |iigo| / 

p and has twice a smaller amplitude. Larger deviations 
from 2(7o lead to progressively decreasing shifts. 

These results are illustrated in FiglSj and show a 
good correspondence to the numerically determined QB 
stability. Summing up, the spatial Foirier components 
p e [2qo — 2, 2qo + 2] of {^n} are decisive for the go-QB 
stability. The dependence is notably different and much 
more complicated than the possible "naive" expectation 
that harmonic inhomogeneities with p ~ qa will most ef- 
fectively stabilize or destabilize go-QBs. A remarkable 
fact is the sensitive dependence on the phase of the im- 
purity harmonics: even for a fixed p opposite shifts in 
the threshold occur. Presumably, this is the consequence 
of the deformation of the linear spectrum due to inho- 
mogeneities, as AVq„,q„ is, actually, the difference in the 
frequency shifts of linear modes In its turn, this is 
determined by the boundary conditions, which also affect 
the nonlinearity induced interaction. It clearly highlights 
one of the future directions of study. 

These findings suggest a possibility of controlling the 



0.265 



0.26- 



0.255 



0.25 



0.245 - 



0.24 




10 



15 



20 



25 



30 



FIG. 3: The instability threshold for QBs with go ^10,N = 
32 and k„ = i cos ^ '^Pt^+y^) . Dash-dotted and marked 

lines are numerical results for solid lines are analytical 
estimates 



energy flow between modes. Indeed, by imposing a 
proper spatially periodic modulation of the linear cou- 
pling one can destabilize certain QB exictations and (i) 
speed up equipartition or (ii) stabilize others, where the 
energy will be radiated to and trapped. New QBs may 
also be destabiHzed to arrange a further energy flow. Ex- 
perimentally, in miniature mechanical lattices inhomo- 
geneities can be created, for example, by laser heating, 
either as harmonic or spot impurities, like it was designed 
to control discrete breathers relocation in cantilever ar- 
rays In optical lattices one can implement the same 
technique that has been recently used for generating dis- 
ordered potentials in studies of the Anderson locaHzation 
of light [3| and matter (BEC) Q waves. 

In summary, we have generaHzed the concept of QBs 
to the case of non-zero disorder and analyzed these non- 
linear vibrational modes in weakly disordered DNLS ar- 
rays. We demonstrated, that QBs remain exponentially 
localized in the mode space and stable, if the disorder 
is sufficiently small. Their stability depends sensitively 
on a particular realization of disorder, and may be en- 
hanced or undermined. The prevaiHng contribution to 
the stability is made by the spatial harmonics of disor- 
der which wave numbers are close to twice of that of the 
QB seed mode. Thus, inhomogeneities design appears to 
be a promising technique of controUing the energy fiow 
between nonlinear modes. We expect these ideas and 
methods to apply to a variety of nonlinear weakly disor- 
dered lattices - and have already demonstrated it for the 
FPU chain (to be reported elsewhere) - to name the dis- 
order coming from other sources (masses, nonlinearities) , 
higher dimensions, and quantum arrays. We believe that 
the results on the nonlinear modes sustainability, stabil- 
ity, and controlling will be widely adopted in experiments 
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and applications. 

We thank S. Flach for stimulating discussions. 
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